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DEFINITION OF SYMBOLS 


Symbol 

Definition 


A 

Propellant tank wall surface area 

Ft 2 

a d 

Distributor area 

CN 

+-> 

hi 

J 

Dimensional constant 

lb£ Ft/Btu 

M w 

Pressurant molecular weight 

lbj^/lb NIol 

T a 

Ambient temperature 

°R 

t l 

Propellant temperature 

°R 

T m 

Ullage mean temperature at cutoff 

°R 

T o 

Pressurant inlet temperature 

°R 

V 

Propellant tank volume 

Ft 3 

Vi 

Initial ullage volume 

Ft 3 

V L 

Propellant volumetric drain rate 

Ft 3 /sec 

d w 

Wall thickness 

Ft 

C P 

Pressurant specific heat 

Btu/Ib m °R 

^pw 

Wall specific heat 

Btu/lb m ° R 

§c 

Dimensional constant 

lb m Ft 


lbf sec 2 

h a 

Ambient heat transfer coefficient 

Btu/ sec Ft 2 0 R 

k 

Pressurant thermal conductivity 

Btu Ft/ sec Ft 2 

P 

Ullage pressure 

lb f /Ft 2 

r 

Propellant tank characteristic radius 
(maximum radius for cylindrical tanks) 

Ft 


IV 



Symb ol 

DEFINITION OF SYMBOLS (CONCLUDED) 
Definition 


R 

Universal Gas Constant 

lbf Ft/ ° 

a 

Gas Compressibility 

( ~ ) 

M- 

Pressurant viscosity, M/L 

lb m /Ft 

Pw 

Wall density, M/L 3 

lb m /Ft 

e 

Time of pressurization 

sec 


R lb Mol 


sec 


v 




PREDICTION OF PROPELLANT TANK PRESSURIZATION 

REQUIREMENTS BY DIMENSIONAL ANALYSIS 

By J. F. Thompson* and M. E. Nein 
George C. Marshall Space Flight Center 


SUMMARY 

an !W1 

Pressurant gas requirements for launch and space vehicles may be 
predicted by analytical models of the pressurization process. How- 
ever, preliminary design studies require a fast and reasonably 
accurate method of predicting without resorting to computer programs. 
Therefore, dimensional analysis of a large number of pressurization 



tanks with an accuracy of ± 10 per cent. 


INTRODUCTION 


Although the most accurate method of predicting pressurant re- 
quirements is with a computer program that has been varified by 
experiments, it is advantageous to have a fast, reasonably accurate 
method to determine the total mass of pressurant gas required with- 
out resorting to the computer. This type of analysis is necessary in 
comparison and optimization studies for preliminary design where 
the number of possibilities to be considered precludes a detailed 
computer analysis of each case. This report presents a single, 
general expression for the total required mass of pressurant; the 
expression was developed by dimensional analysis of the results of 
about 30 pressurization tests and 120 simulations on an IBM 7094 
Computer. 


* Assistant Professor at Mississippi State University Aeronautical 
Engineering Department; formerly with MSFC. 




DIMENSIONAL ANALYSIS OF TEST AND COMPUTER RESULTS 


The total mass of pressurant gas required is a function of the 
ullage mean temperature at cutoff derived with the gas equation of 
state, 


PV M w 

w = 

Total a RT m 


( 1 ) 


Therefore, the total pressurant mass required may be calculated if 
the ullage mean temperature at cutoff can be determined. In the 
most general case, the ullage mean temperature at cutoff will be a 
function of 1 Z system design variables, seven physical properties, the 
mechanical equivalent of heat, and the gravitational constant: 


f(J- g c , 


k, p, c p , 

r, T q , Tj^ , 0^ , V, A, 


a ’ C Pw ’ 


d , p, V. 

W ‘ 1 

* A * a d ^ 

(2) 


The 19 variables can be expressed in six fundamental dimensions, 
length (L), mass (M), time (0), temperature (T), heat (H), and force 
(F). The two dimensional constants, J and g c , are included because 
heat and force can be expressed in terms of the other four dimensions 
The dimensions of each variable are given in Table IV. 

Since any equation representing physical phenomena must be 
dimensionally homogeneous, it must be possible to write equation (Z) 
in a nondimens ional form. Therefore, using tt as a symbol for a 
dimensionless group, equation (Z) may be written as follows: 


*1 = ffa . "3 , TT ) 


( 3 ) 


where the group ttj contains the ullage mean temperature at cutoff. 
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According to Buckingham's theorem, the number of dimensionless 
groups in equation (3) is given by: 


i = n-m 


(4) 


where n = number of variables 

m = maximum number of these variables which will not form 
a dimensionless group 

The quantity m is often equal to the number of fundamental 
dimensions. The total number of variables in equation (2), including 
J and g c , is 22, and the number of fundamental dimensions is six. 

By trial and error, it was determined that the maximum number of 
variables in this case that will not form a dimensionless group is 
indeed six, the number of fundamental dimensions. Therefore, 
according to equation (4) there will be 16 dimensionless groups in 
equation (3). 


Dimensionless Groups 

Because the maximum number of variables that will not form a 
dimensionless group is six, the 16 dimensionless groups may be 
determined by choosing six variables to be common to all groups and, 
in turn, adding each of the remaining 16 variables to the first six to 
form 16 dimensionless groups. Although any six variables could be 
chosen as the six to be common to all groups, intuition leads to the 
choice of the two dimensional constants (J and g c ), the pressurant 
molecular weight (M ), thermal conductivity (k), and viscosity (p) 
and the propellant tank characteristic radius (r). Adding the ullage 
mean temperature at cutoff (T ) to the six common variables yields 
the first group: 


a b c d e f 
Tr 1 =Jg Mk p r 
c w 


m 


(5) 
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now, substituting the dimensions: 




"FL" 

a 

'ML 

D 

c 

M 

H 1 

d 

M 

e 

L 

1 

'1 


. H J 


Lf 6 2 j 


- L0T - 


LL0 - 



L T J 


g 


( 6 ) 


Since the quantity on the left of equation (6) is dimensionless, the 
exponent of each of the six dimensions on the right must be zero. 

This condition determines six simultaneous equations, 


L L J 

0 

= a + b - d 

M 

0 

= b + c + e 

0 

0 

= -2b-d-e 


0 

= -d + g 

b] 

0 

= »a + d 


1 0 

= a - b 


(7) 


Simultaneous solution of these equations yields the following: 


a ~g 

bJD 

II 

"d 

b~g 

e- ”3g 

c = 2g 

£=-4g 


4 



Because the value of g is arbitrary, it can be taken as unity and 
yields the following: 


a = 1 d = 1 

b = 1 e = -3 

c = Z f = -4 


Substituting these values into equation (5) yields the first group : 




g c M w k 


L 3 r 4 


T 

m 


( 8 ) 


The second group tt 2 is found in the same manner by adding the 
pressurant specific heat C p to the six common variables. Each of 
the remaining variables is added one at a time to the six common 
variables, and the remaining dimensionless groups are determined 
in the same manner as the first two. Then, 


"2 


Cp p 
k 


(9) 


TT3 


gc 


M 2 k T 
w o 

p 3 r 4 


tt 4 


gc 


M 2 k 
w 

|jl 3 r* 



( 10 ) 


( 11 ) 


^5 


p r 


M 


w 


( 12 ) 


"6 = 


V 


(13) 


5 



*7 


*8 


tt 9 


*10 


*11 


*12 


*13 


*14 


*15 = 


*16 


Cp W H; 


p 3 




iw_ 

r 


Vi 


gc p 

p 2 r 


h a r 
k 



M 2 k' T 
w a 


M w V|. 

n r 4 


A P 

r 2 


(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23 
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Now, the groups ttj , Tr 4 , and tt i4 may be divided by tt 3 with no 
loss of generality, because none of the 16 groups are eliminated in 
the process. Thus Hi , tt 4 , and tt 14 maybe replaced with tt* , tt 4 , and 
ttj 4 , respectively, where 


/ 

*1 


IT 4 = 


tt'i 


14 


7T = TT 7T 
11 11 1 



_ T m 

■n-3 

To 

tt 4 

_ t l 


T o 


_ T a 

^3 

To 

P *11 

may be replaced 


Kc. MwP Vi 

Z 

p 2 r 4 


( 24 ) 


This change occurs because the initial ullage mean temperature is 
proportional to the product pV^ . 

It is possible to reduce the number of dimensionless groups 
necessary in this case by realizing that the wall specific heat, density, 
and thickness can enter the problem only in the combination Cp W p w d w - 
which is the wall heat capacity per unit area. Therefore, the groups 
7t 8 , tt 9 , and Tr 10 may be combined into a single group *rr 8 . 


/ 

TTg - TT 8 TT 9 TT i0 


M w k 

p r 2 (C p d ) 
^ ' pw r w w 7 


( 25 ) 
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Similarly, it is possible to combine ir 15 and tr l6 . 


Thus, 


*15 = 


*15 

*16 


ft , 

A£)P r 2 


( 26 ) 


Further, the propellant tank volume and wall surface area should 
enter the problem as the ar ea -to-volume ratio (A/V). Therefore, the 
group it 6 and tt 7 may be combined into a single group . 


/ 

*6 - 


*6 


Ar 

V 


(27) 


But since the approximate area-to-volume ratio is given by 


A ^ 2tt r 1 _ 2_ 

V~ = it r 2 1 “ r (28) 


the group is nearly a constant and may be eliminated. 

Finally, since the Prandtl number (tt 2 ) does not vary greatly, this 
group may also be dropped. Thus, the number of dimensionless groups 
necessary in equation (3) was reduced by six. 

Since the propellant temperature is the lower limit of the ullage 
mean temperature, it is logical to replace all temperatures with 
temperature differences above the liquid propellant temperature (T t ). 

Equation (3) may now be expressed in terms of the following ten 
dimensionless groups (the group numbers having been changed to run 
from one to ten): 


l 



(29) 
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^z = 


tt 3 = 


tt 4 = 


tt 5 = 


^6 = 


tt 7 - 


TTg = 


tr 9 = 


TTlO 


J g C M w k (Tq 


Jjj 


(i 3 r 4 


- T t 


T, 


L 


(u~ Op 


M 


w 


Hyfc 

I* ri < C p'» P w d w> 


g r M w p Vj 


( 30 ) 


( 31 ) 


( 32 ) 


( 33 ) 


( 34 ) 


&c Ma P (35) 

fj . 2 r 

h a r (36) 

k 

T a - TL- ( 37 ) 

T o - T L 

Mw Vt. ( 38 ) 

Aq p. r 2 
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Curve Fit of Dimensionless Equation 

These ten dimensionless groups [equations (29) through (38)] 
can be used to correlate the results of tests and computer runs 
according to equation (3). In most cases equation (3) would be written 
in the following form: 


_ 07 6 « i X T o- i 

7 h - O' *2 tt 3 tt 4 ir 5 tt 6 tt 7 (39) 

However, in this case it is necessary to satisfy certain boundary 
conditions that cannot be satisfied by equation (39). The ullage mean 
temperature at cutoff must remain finite and not equal to zero as the 
ambient heat transfer approaches zero. This boundary condition cannot 
be satisfied by equation (39), but it can be satisfied if the functional 
dependence on ng and tt 9 is exponential. Also, as the distributor 
Reynolds number tt 10 approaches zero, the heat transfer in the tank 
approaches free convection. Therefore, the boundary condition of 
finite, non-zero mean temperature, when ir 10 is zero is imposed, 
dictates an exponential functional dependence on tt 10 . Thus, these 
boundary conditions can be satisfied by writing equation (3) in the form 


"i * ffi W 3 4 4 ttj 4 4 M (40) 

The coefficients and exponents in equation (40) were evaluated by a 
curve fit to the data from the computer runs and tests. It was found 
that the data could be correlated by this equation if the coefficients 
ni and #3 , in the exponentials were taken as functions of tt 2 and 1T3 ; 


/ 00 4» 

a 2 - ot 2 TT2 ^3 


(41) 


t 

a 3 


4> 

a 3 "2 


(42) 
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Equation (40) then becomes 


3 y 6 e i X 

TTi -QTi TT 2 TT 3 TT4 TT 5 TT 6 TT7 



GO 

^2 


TT3 


T 

TTg 



( 43 ) 


where all coefficients and exponents are constants. 

It was convenient to divide several groups of equation (43) by powers 
of ten to obtain numbers more easily handled. Therefore, the final 
equation used in the curve-fit was 



From the curve-fit, the following values were obtained for the 
coefficients and exponents in equation 44: 


ai = 0. 4Z4 

a z = 0 . 00210 

= -0. 0292 
3 = -0. 1322 

y = -0. 1688 

5 = -0. 1146 

e = 0. 0780 


£ = 0. 01416 
X = 0. 0620 
u = 0.415 
vj> = 1. 174 
T = 0. 765 
<{> = 0. 1510 
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This equation is general and is capable of predicting the ullage mean 
temperature, and thus pressurant mass at cutoff within 10% for 
cylindrical tanks with rounded bulkheads* FIG 1 shows total pressurant 
requirements obtained by various investigators for a wide range of tank 
sizes and system parameters compared with the pressurant weights 
calculated by equation (44). Excellent agreement is obtained over the 
entire range of conditions for hydrogen and oxygen pressurization. 
However, the equation is limited in its application to conditions of 
constant ullage pressure, pressurant inlet temperature, and ambient 
heat transfer. The studies indicated that the equation is inaccurate 
at inlet temperatures less than 100°R above the saturation temperature, 
at ullage pressure below propellant saturation, and for very short 
expulsion times of less than 50 seconds. The restriction to cylindrical 
tanks can be removed by proper choice of the characteristic tank radius. 

Studies have shown that the characteristic tank radius for oblate 
spheroids, used in equation (44), should be about two-thirds of the 
maximum tank radius. This assumption is theoretically justified, 
because a cylinder having the same volume and surface area as an 
oblate spheroid has a radius equal to 0. 63 times its maximum radius. 
Further test data and analytical studies are necessary to select the 
characteristic radius for other geometries. 


1Z 



Due to the dimensionless nature of equation (44), it is not restricted 
to any particular propellant, pressurant, or tank size as indicated in 
FIG 1. Although an uninsulated propellant tank was assumed in the 
development of this equation, FIG 1 shows good agreement with test 
results obtained with vacuum jacketed liquid hydrogen tanks. 

To simplify the use of equation (44), the pressurant and propellant 
properties for the case of liquid oxygen pressurized by oxygen and 
liquid oxygen pressurized by helium were substituted in equation (44), 
and the following equations were obtained. Since these equations are 
dimensional, they are applicable only to the case indicated. For liquid 
oxygen pressurized by oxygen, 


T m - 164 

T q - 164 


3. 33 (T q - 164)“°* 297 


0.1395 „ 0* 01416 

r V: 


^Pw d w) 


-0. 0780 0. 0762 -0. 1146 

P 6t 


exp 


0. 001568 (T 0 - 164) r 


-0.304 -0. 895, 0. 765 


ha (T„ - 164) 


* exp 


-120. 9 (T q - 164) 


-0. 574 -2 


. 604 / Vl\ 

Ud/ 


(45) 


and for liquid oxygen pressurized by helium, 


T m - 164 


To - 164 


= 3. 10 (T 0 - 164) 


-0. 304 0. 1395 Tr 0. 01416 


r ^ 


, -0. 078 0. 0762 -0. 1146 

‘( C n P 0 T 


w 


exp 


0.0000681 (T - 164) 


-0. 1650 -0. 895 , 0. 765 


ha (ta -164) 


-5.04 (T 0 -l64)-°--r- 2 -^ ^ 


(46) 


. exp 
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The units of all variables in equations (45) and (46) must be those given 
in the Definition of Symbols. For other combinations of propellant and 
pressurant, equation (44) must be used. After the ullage mean tempera- 
ture at cutoff is calculated, and using equation (44), (45), or (46), if 
applicable, the total mass of required pressurant gas may be calculated 
from equation (1). 

In designing a launch or space vehicle pressurization system, vehicle 
parameters such as tank volume, engine flowrate, tank material, etc. , 
determined by vehicle mission profile, are fixed input values. How- 
ever, there are various controllable parameters in a pressurization 
system that can be used to optimize the system without affecting 
basic vehicle characteristics. The relative significance of various 
parameters on pressurant requirements has, therefore, been investi- 
gated under another study program. The results of these studies 
excerpted from reference 1 are presented in FIG 2. From a central 
origin, representing a reference condition (SATURN V, S-IC Stage) for 
all parameters, the increase (+Y) and decrease (-Y), of the ullage 
mean temperature at cutoff is shown as a function of variation of the 
parameters on the abscissa. The parameters were varied over a 
range expected for vehicle design. Thus, pressurant inlet temperature 
can increase or decrease by a factor of 2 from the reference condition, 
pressure by a factor of 3, tank radius by a factor of 2, expulsion time 
by a factor 3, etc. It was indicated that the pressurant inlet temperature 
exerts the greatest influence on the ullage mean temperature. Diminish- 
ing return of this effect did not exist within the range of investigation 
(530°R to 1200°R). The mean temperature increased as the ullage 
pressure was increased and also as the tank radius was increased. 
Increasing the tank wall thickness, heat capacity, or density caused a 
decrease in the mean temperature. The pressurant distributor flow 
area (A-q) that controls the gas-to-wall forced convection heat transfer 
coefficient had a significant effect on the mean temperature when Ajj 
was reduced, but no effect at all when flow area was increased. This 
indicates that the pressurant inlet velocity for the reference systems 
was chosen at an optimum point. FIG 2 also indicates that helium 
pressurant must be introduced into a tank at a temperature 1. 1 times 
higher than oxygen pressurant to obtain the same ullage mean tem- 
perature. 
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CONC LUSIONS 


An equation derived by dimensional analysis provides a reasonably 
accurate method for prediction of pressurant requirements for 
cylindrical LOX and hydrogen propellant containers. This method 
is advantageous for preliminary design and optimization studies where 
the use of large computer programs becomes excessive in cost and 
time. 
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